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Abstract 

The results reported in this paper describe some of the main 
flow characteristics and NOx production results which 
develop in the mixing process in a constant cross-sectional 
cylindrical duct. A 3-dimensional numerical model has been 
used to predict the mixing flow field and NOx characteristics 
in a mixing section of an RQL combustor. Eighteen configu- 
rations have been analyzed in a circular geometry in a fully 
reacting environment simulating the operating condition of 
an actual RQL gas turbine combustion liner. The evaluation 
matrix was constructed by varying three parameters: 1) 
jet-to-mainstream momentum-flux ratio (J), 2) orifice shape 
or orifice aspect ratio, and 3) slot slant angle. The results 
indicate that the mixing flow field and NOx production 
significantly vary with the value of the jet penetration and 
subsequently, slanting elongated slots generally improve the 
NOx production at high J conditions. Round orifices produce 
low NOx at low J due to the strong jet penetration. The NOx 
production trends do not correlate with the mixing 
non-uniformity parameters described herein. 

Nomenclature 


Ajk cross-sectional area of the control volume at i and j 
location, m 2 

An, duct cross-sectional area, also A^, m 2 
ACd effective orifice area, m 2 

Ar area ratio (jet/mainstream)= ACd/ An, 

AMIX area weighted deviation from equilibrium, Eq-3 
B area determined half width of the distribution func 
tion 
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c 

constant of proportionality in Eq-6 

d 

diameter of the orifice 

DR 

density ratio (jet/mainstream) 

DP/P 

total pressure loss across the mixing wall, % 

f 

non-dimensional equivalence ratio, Eq-2 

f/a 

fuel to air ratio 

J 

momentum-flux ratio (jet/mainstream) = M 2 /DR, 
also (MR) 2 /[(DR)(ACd/A m ) 2 ] 

m tot 

overall mixer mass flowrate, kg/sec = nv, + 

M 

mass-flux ratio (jet/mainstream) = DR VWU mam 
mass flow weighted deviation from equil, Eq-4 

MMDC 

MR 

mass flowrate ratio (jet/mainstream) 

n 

optimum number of orifices / row, Eq-6 

P 

total pressure, atm. 

R 

radius of the mixing section, m 

r 

radial distance from the centerline of the mixer, m 

S 

spacing between orifice centers at the circumference 
of the can 

T 

temperature, K 

u 

local axial velocity, m/sec 

iimain 

approach mainstream axial velocity, m/sec 

VR 

velocity ratio (jet/mainstream) =Vj et / U ma in 

Vje, 

radial velocity of the jet, m/sec 

Z 

arc distance in the tangential direction evaluated at 
r=R; z equals 0 at midplane between orifices m 

X 

axial distance from the leading edge of the orifice, m 

5 

value of the cumul volume fraction at f^uu, Fig-4 

P 

fluid density, kg/m 3 

9 

equivalence ratio (f/ajjocai / (f/a) sto i 

Subscripts 

equil 

equilibrium 

i 

axial vector direction 

j 

radial vector direction 

k 

tangential vector direction 

m 

mainstream, also (main) 

RZ 

rich-zone 

stoi 

stoichiometric 
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Introduction 

In recent years, the design and development of gas turbine 
engines for the aeropropulsion and ground based power 
generation systems has been channeled towards decreasing 
the gaseous emissions without adversely affecting the system 
performance. The environmental effects of both carbon 
monoxide (CO) and oxides of nitrogen (NOx) have been 
investigated extensively for many years and their potential 
hazards are increasingly becoming a more world wide sensi- 
tive issue. The trend in gas turbine engine design advance- 
ment has been channeled towards increasing both engine 
pressure ratio and rotor inlet temperature levels in order to 
increase the overall thermodynamic cycle efficiency. How- 
ever, this evolution has been adversely affecting both engine 
durability and gaseous emissions, especially that of NOx, 
since the production of NOx is generally a function of the 
combustion system design, hot section temperature distribu- 
tion, residence time, and localized fuel/air mixture. 

The development of an efficient and low-emission combus- 
tion system demands increased insight into combustion 
chemical kinetics, efficient air/fuel mixture, and advanced 
materials and cooling techniques. Current gas turbine com- 
bustion liner technology employs a single-stage combustion 
process in which both fuel and air are admitted into a con- 
trolled mixing chamber and are allowed to react. However, 
with the advent of higher temperature operating conditions, 
the reduction of NOx becomes a very difficult task to accom- 
plish using single axial staged combustion. Therefore, 
alternative combustion methods are being explored by the 
primary gas turbine engine manufacturers and other research 
organizations. Two of the main low NOx designs being 
developed are rich burn / quick mix / lean bum (RQL), and 
lean premixed prevaporized (LPP). Both have advantages 
and disadvantages regarding both operation range and hard- 
ware complexity. 

This paper focuses on the performance of the mixing section 
of an RQL combustor. The successful performance of this 
combustor depends on a quick and efficient mixing of the 
rich zone combustion products with compressor discharge air 
to effectively reduce the overall equivalence ratio from about 
1.8 to about 0.5. This process must be accomplished with a 
minimum transient time (near equivalence ratio of 1) where 
most of the NOx is produced due to the high resulting tem- 
perature levels and oxygen availability. 

An experimental and analytical effort (e.g. Bain, Smith, and 
Holdeman, 1992, 1993, and 1994; Doerr and Hennecke, 
1993; Hatch, et al„ 1992a and 1992b; Howe, et al., 1991; 
Kroll. et al., 1993; Liscinsky et al.. 1992; Liscinsky, True, 
and Holdeman, 1993 and 1994; Oechsle, Mongia, and Holde- 
man, 1992. 1993, and 1994; Smith, Talpallikar, and Holde- 
man, 1991; Sowa, et al., 1994); Talpallikar, et al., 1991; 
Vranos, et al., 1991; and Zhu and Lai, 1992) is underway to 
study and identify the critical design and flow parameters 
affecting the mixing effectiveness. 


In this study, a 3-D numerical tool is used to predict the 
performance of an RQL mixing section. The 3-D flow 
structure of the mixer can be modeled in detail and predic- 
tions can be obtained with a host of scalar and vector quanti- 
ties to accurately evaluate the mixing flowfield. In addition, 
a relatively large number of configurations can be easily 
analyzed to predict the trends of both mixing uniformity and 
NOx production in a generalized jet in a crossflow problem. 
In this study, the flowfield of several different mixer configu- 
rations are evaluated and compared with their corresponding 
NOx predictions. For this purpose, 18 different mixer con- 
figurations were analyzed with the 3-D numerical tool and 
the localized NOx and overall mixer NOx production flow- 
rates were calculated from each case. The mixer configura- 
tions include the following orifice shapes: a) round holes, and 
b) elongated slots with aspect ratios of L/W=4 and 8 and 
different slant angles from 0 degrees (aligned with the flow) 
to 90 degrees (transversely oriented to the flow). The jet to 
mainstream momentum-flux ratio was varied from 25 to 80. 
An analysis was further carried out to evaluate and quantify 
the mixing flowfields in order to rank mixing and NOx 
production configurations. 

Mathematical Model 

3-D Flow Model 

A production 3-D combustor code, COM-3D (Bruce, Mon- 
gia, and Reynolds, 1979) is used that solves the turbulent 
reacting flow transport equations using the SIMPLE algo- 
rithm of Patankar and Spalding (Patankar, 1980). This 
program simulates turbulence by the two-equation k-e model 
(Launder and Spalding, 1974), and combustion following 
vaporization is determined by a four-step chemical reaction 
model based on Arrhenius and modified eddy breakup con- 
cepts. The transport equations for all dependent variables are 
of the following form as shown in Eq-1 

div[p r u % - (Peff/Pr) grad(^)] = (1) 

where p r is the mixture density, u is the velocity, mff is the 
effective turbulent viscosity, P r is the effective 
Prandtl/Schmidt number, and S^ is the source term for the 
variable £. The following variables are computed by 
COM-3D; 1) axial, radial, and swirl velocity components; 2) 
specific enthalpy and temperature; 3) turbulence kinetic 
energy and dissipation rate; 4) unbumed fuel, CO, H 2 , inter- 
mediate fuel, and composite fuel mass fractions; and 5) fuel 
spray trajectory and evaporation rate. 

The computational effort is significantly reduced by model- 
ing a sector of the mixing section comprising a single orifice. 
Therefore, the shape of the sector was dependent on the 
number of orifices equally spaced in the circumferential 
direction. It should be noted that all configurations shown in 
this paper contain 8 orifices/row and each yield a 45 degree 
sector computational domain. Periodic boundary conditions 
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were applied in the circumferential direction. No-slip and 
adiabatic boundary conditions were applied at the outer wall 
defining the inside wall of the mixing section. Zero-gradient 
boundary conditions were applied at the center axis. Axial 
gradients at the exit boundary condition were assumed zero. 

NOx Model 

The NOx model described herein was developed by Rizk and 
Mongia (1993). Because the NOx formation in the combus- 
tion liner is significantly affected by the details of the front 
end of the RQL and the subsequent admittance of air into the 
various downstream zones, the combustion liner needs to be 
divided into a number of regions for modeling purposes. The 
hybrid modeling technique therefore consists of using the 
3-D simulation results obtained with COM-3D such as gas 
flowrate, flow averaged temperature, fuel/air ratio, in addi- 
tion to the turbulence characteristics to accurately describe 
the flow nuances affecting the NOx production. 

The input flow field exiting the rich zone and entering the 
mixing zone is assumed to be fully reacted and in equilib- 
rium. Species mole fractions for (CO, CO 2 , H 2 0, and H 2 ) 
were determined based on the given rich-zone equivalence 
ratio at chemical equilibrium at the prescribed operating 
conditions using JP-5 for typical fuel properties. This 
assumption was necessary in order to duplicate the inlet 
flowfield to all the mixing configurations analyzed herein. 
The four specified species used in the inlet mainstream 
condition correspond to the four-step chemical reaction used 
in COM-3D. 

The hybrid model consists in grouping the results obtained 
from a typical COM-3D grid with about 50,000 to 80,000 
nodes into about 2,000 larger sub-volumes comprised of 
several computational nodes representing the entire mixing 
section. The reacting flow model results can therefore be 
accurately represented and modeled using significantly less 
computational effort. For the NOx predictions only, the 
model of the mixing zone extends to x/R=5. This was imple- 
mented to prevent any unrealistic discontinuities in the 
predictions near the real interface between the mixing and 
lean zones (x/R=l) where x is the downstream axial distance 
from the leading edge of the orifice and R is the mixer radius 
evaluated at the outer wall. 

Geometric Configuration 

In this study, the mixing section was modeled as a constant 
diameter cylindrical duct with a single row of equally spaced 
orifices. The outer wall diameter is 3 inches (0.076 m) and 
the axial length of the mixing section extended from 
x/R=-1.4 to x/R=6. Sufficient axial distance was provided 
both upstream and downstream of the orifice to minimize the 
impact of the inlet and exit boundary conditions on the 
calculated flow structure in the primary domain of interest 
which is (0<x/R<l). The downstream limit (x/R=l) of the 


mixing zone has been defined arbitrarily by this and many 
other authors listed in the reference. 

The computational grid domain was typically discretized into 
50,000 to 80.000 finite control volumes generally arranged 
with about 70 nodes in the axial direction, and 30 nodes in 
both the radial and tangential directions. The grid was 
typically denser near the orifice and near the outer wall to 
resolve the high velocity and temperature gradients resulting 
from the inlet of the crossflow jet. An orthogonal view of a 
typical grid arrangement is shown in Figure- 1. The grid is 
normally configured to allow smooth progressive volume 
change between adjacent control volumes to help speed up 
the convergence of the solution. In the 3-D numerical model, 
a secondary grid (staggered grid) is interpolated from that 
shown in Figure- 1 to obtain the boundaries for the vector 
quantities. The staggered grid becomes the control surfaces 
of the micro control volumes that constitute the inner volume 
of the mixing section. 

The geometric configuration of the mixer showing a slanted 
slot is also shown in Figure- 1. A total of 18 circular and 
slanted slot hole configurations were analyzed as shown in 
Tables- 1 and 2 and are also shown graphically in Figure-2. 
The blockage is defined as the circumferential projection of 
the orifice divided by the spacing between the orifice centers, 
and the T.E. x/R is the non-dimensional axial location (x/R) 
of the trailing edge of the orifice with respect to its leading 
edge (x/R=0). Both these parameters are shown in Tables- 1 
and 2. The matrix of slanted slot orifices summarized in this 
report include shape variations with an aspect ratio (L/W = 
long/short) of 4 and 8, and slant angles of 0, 22.5, 45, 67.5, 
and 90 degrees with respect to the mainstream flow direction. 
Note that the nomenclature that identifies the configurations 
range from 19 through 36. Configurations 1 through 18 have 
previously been reported by Oechsle, Mongia, and Holdeman 
(1993), and they are the non-reacting counterparts of the 
same geometric configurations. 

In this study, the following main parameters were kept 
constant throughout this analysis. The control of these 
parameters is essential for a valid comparison of both mixing 
non-uniformity and NOx production between the appropriate 
configurations. 

1) rich zone equivalence ratio <{> ra = 1.80 

2) lean zone equivalence ratio = 0.416 

3) overall mixer pressure = 14.1 atm. 

4) jet temperature = 950 K 

5) mainstream temperature = 2164 K, from the chemical 

equilibrium code (CEC) developed by NASA 

6) mixer diameter = 0.076 m 

7) jet to mainstream mass flow rate MR = 2.96 

8) jet to mainstream density ratio DR = 2.28 

9) number of equally spaced orifices per row = 8 

The following were allowed to vary: 
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1) jet to mainstream momentum-flux ratio J from nominal 
values of 25 to 80 

2) overall mixer total pressure drop between 1.3 to 5% 

3) orifice ACd depending on the value of J in order to 
maintain the same MR 

4) orifice blockage and T.E. x/R values depending on the 
orifice shape and orientation 

5) jet velocity depending on the value of J 

Modeling Specifications 

The species mole fractions, temperature, and velocity profiles 
exiting the rich zone and entering the mixer were assumed to 
be uniform across the inlet cross-section of the mixing sec- 
tion. The air jet flow was characterized by a radial, uniform 
flow across the orifice effective area. COM-3D is not a body 
conforming code therefore accurate modeling of the orifice 
shape was obtained by defining the orifice with 80 to 150 
control surfaces. In this way the stair-stepping approxima- 
tion in the slanted and round contours could be maximized. 
The assumption of uniform mass injection/area for the orifice 
air entry is applied in the mathematical model in all the 
analyzed configurations. The turbulence kinetic energy of 
the mainstream and jet flows were 0.3% of the square of the 
mean velocities. The turbulence length scales of the main- 
stream flow were 2% of the mixer diameter, and the turbu- 
lence length scale of the jet was of the order of the orifice 
diameter. The results from COM-3D were post processed 
into about 2000 sub-volumes typically 20 in the axial direc- 
tion, 10 in the radial, and 10 in the circumferential direction. 
This input was used to generate the NOx predictions based on 
the flow characteristics. 

A typical numerical solution took about 500 iterations for full 
reacting flow convergence with overall mass flow residuals 
of 0.05% of the total mixing section mass flowrate. All 
solutions were obtained using the Cray C-90 and a typically 
converged solution took about 2 hrs of CPU time. 

Analysis 

The mixing performance for all configurations analyzed in 
this study were ultimately evaluated at x/R =1. It is however 
recognized that the performance throughout the mixing 
section volume of interest (0 < x/R < 1) should also be con- 
sidered since complex structures are present in the flow field 
especially near the entry of the jet. Two different methods 
were used to analyze the numerical results: observational 
analysis and statistical analysis. Both of these are described 
in detail below. 

Observational analysis 

Normalized equivalence ratio contour distributions, and 
velocity vector plots are represented for each configuration. 
The equivalence ratio results are presented as normalized 
values with respect to the overall differential between the 
mainstream flow and the jet flow inlet values. The param- 


eter (0 is defined in Eq-2 . Note that the parameter (f) is a 
conserved scalar. 

f jk~ i et ( 2 ) 

(j) main <j) jet 

The value of (f) varies from 0 to i, where 0 is the equivalence 
ratio value of the unmixed jet and 1 is the equivalence ratio 
value of the mainstream flow. Note that f = 1-9, where 0 is as 
defined previously (Holdeman, 1993) and used elsewhere 
also. The equilibrium (0 value was calculated in a purely 
adiabatic system at any location downstream of the jet injec- 
tion. The equivalence ratio and velocity plots as shown in 
the example plot in Figure-3a are given in the axial-radial 
plane through the center of the orifice where the axial direc- 
tion is in the same direction as the mainstream flow direction 
(flow towards the turbine). Both axial and radial directions 
are non-dimensionalized with respect to the mixer radius. 
The x/R=0 location denotes the leading edge of the orifice 
and a r/R location of 1 is the outer wall location. The view is 
such that the left side of the plot is the upstream side of the 
mixer. 

Three-dimensional orthogonal views of the local NOx pro- 
duction in gr of NOx / {sec m 3 }, fuel/air ratio, and tempera- 
ture are shown in this paper (see example plot in Figure-3b 
for the NOx production only). The 3-D orthogonal view 
plots show four different planes that essentially lengthwise 
split the mixer domain in half. The planes shown are: a) the 
radial-tangential planes at slightly upstream of x/R=0, b) 
slightly downstream of x/R=l, c) axial-radial plane through 
the center of the orifice, and d) axial-tangential plane at the 
outer wall of the mixing section. In order to eliminate any 
3-dimensional ambiguity in the interpretation of these 3-D 
plots, it is necessary to understand that the viewer is looking 
up at the inside surface of the outer wall of the mixing sec- 
tion and half of the orifice can be seen since the axial-radial 
plane slices the orifice. It should also be noted that the 
viewer is located downstream of the mixing section. In these 
plots, the flow moves from left to right. 

The interpolated iso-surface of lCf 3 gr NOx / {sec m 3 } is also 
shown in a different perspective as compared to the previous 
plots (see example plot in Figure-3c). These 3-D orthogonal 
plots show the mixing section (0<x/R<l) looking down- 
stream from the upstream end of the mixer. This edge closest 
to the observer is shown to eliminate the ambiguity of the 
3-D perspective. Note that the outer contour (showing the 
pie section of the mixer) is slightly smaller than the actual 
mixer (the outer wall is located at r/R=l) since it is described 
by the center points of each of the outer most sub-volumes in 
the NOx analysis. Figure-3d shows the mixer domain and 
sub-volume breakup used in the NOx calculations. The 
intersecting points correspond to the center of the 
sub-volumes. The true 3-D location of the orifice is also 
shown. The flow is moving from left to right, and the coordi- 
nates have been normalized in the x, y, and z directions. 
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Statistical analysis 


Volume Fraction 


This analysis method involves the detailed characterization 
of the mixing flowfield using several different statistical 
techniques. This technique is especially useful for the 
numerical results since a copious amount of data are obtained 
from the numerical model with a host of different scalar and 
vector quantities at several thousand locations in the flow- 
field. The statistical techniques reported herein therefore 
characterize the mixing non-uniformity only. Three different 
statistical methods are reported in this paper 

a) The performance of the mixing section at the 
radial-tangential planes at x/R=l was quantified by 
using area weighted planar deviation parameter. The 
smallest deviation with respect to the equilibrium value 
indicates the best mixing configuration. This parameter 
(AMIX) is described in Eq-3. It is also important to 
note that AMIX does not correct the mixing 
non-uniformity for the bias introduced in the region of 
the mixer where the air is being injected through the 
orifice. However, AMIX is only evaluated at the exit 
of the mixer x/R=l and is therefore applicable since this 
plane is downstream of the air injection in all 18 con- 
figurations summarized in this report. The 
area-weighted non-uniformity results are shown in 
Table-3. 



b) The mass flow weighted planar deviation parameter 
(also evaluated at x/R=l) as defined in Eq-4 is also 
used to evaluate the mixing region. Note that this 
parameter is similar to the area weighted parameter in 
Eq-3, with the added density and velocity weighting 
terms. The mass flow weighted mixing non-uniformity 
results for all configurations are shown in Table-3. The 
smallest MMIX value with respect to the f^i corre- 
sponds to the best mixing uniformity. 


MMIX = 





0 * ~ 0«1 
0 oiain~ 0 jet 


if 


(4) 


c) The flow field was also evaluated by performing a 
numerical volume integration throughout the mixing 
section of interest (0 < x/R < 1) as shown in Eq-5. 


fhj 


Volume ^ 


bin 


fmain 


' Volume . 


ijk 


(5) 


The volume integration parameter allows the analysis of 
the entire flow field mixture which is more descriptive 
of the overall flow phenomena in the mixing section as 
compared to a planar deviation analysis shown in Eqs-3 
and 4. This volume integration was only performed on 
the equivalence ratio distributions. The entire range of 
the normalized parameter f (from 0 to 1) was 
sub-divided into 200 equal size bins and the volume of 
the computational control volumes corresponding to the 
value of (f) at a certain bin size (fj to f i+Af ) was inte- 
grated as shown in Figure-4. The integrated volume in 
each bin was normalized based on the entire analyzed 
mixer volume, thus obtaining the normalized volume 
fraction. The volume fraction for each bin was plotted 
in the ordinate with the corresponding value of f on the 
abscissa forming a histogram plot. 


The cumulative volume fraction shown in Figure-4 is the 
integrated histogram for the range (0<f£l) and used to 
evaluate the shape of the volume fraction histogram by 
the definition of B(+) and B(-). This evaluation is 
similar to that used by Oechsle, Mongia, and Holdeman, 
(1992). The value of B(+) is the "area determined" 
distribution half width above f^i such that the inte- 
grated area under the histogram above is 1/2 that 
of the overall histogram area above fg^. The same 
applies to B(-) for the area below The definitions 
of B(+) and B(-) are shown pictoriafly in Figure-4. The 
most uniform mixing results yields lowest B(+) and 
B(-), meaning that the histogram distribution width 
effectively collapses towards the equilibrium value of 
(0. The histogram shape value SUM-B (which is the 
value reported in this paper) is defined as the sum of 
B(+) and B(-). 


The summary of both mixing uniformity and NOx production 
flowrate trends for the 18 configurations are described in the 
following sections: 


1 ) effect of orifice shape on mixing and NOx 

2) effect of J on mixing and NOx production, 

3) effect of slot slant angle on mixing and NOx 

4) configuration optimization at different J, and 

5) correlation between NOx production and the mixing 
flowfield 
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Results and Discussion 

Effect of orifice shape on NOx production and mixing 
non-uniformitv 

Note that the corresponding non-reacting distributions appear 
in Oechsle, Mongia, and Holdeman (1992) and correspond- 
ing experimental equivalence ratio plots and NOx inference 
distributions appear in Hatch, et al. (1992a and b). The 
mixing non-uniformity and NOx production were compared 
for the round orifices, 45° slots with L/W=4, and 45° slot 
with L/W=8 for J=25 (configuration numbers 19, 20, and 24 
respectively), at J=52 (configuration numbers 23, 29, and 
25), and J=80 (configuration numbers 22, 30, and 26). These 
configurations are shown in Figure-2. 

J=25 

The results at J=25 shown in Figures-5a and 6 indicate that 
the round orifice produces less NOx for eight orifices at this 
momentum-flux ratio. Figure-5a shows the results of all 
three statistical parameters AMIX, MMIX, and SUM-B 
depicting the mixing non-uniformity. Figure-6 indicates the 
flowfield development throughout the mixing section with 
the normalized equivalence ratio on the left column and 
velocity vectors on the right. As mentioned previously, the 
lowest values are desired for all three parameters. In addition 
to the mixing, the cumulative NOx production in gr/sec up to 
x/R=l is also plotted. The cumulative NOx value at x/R=l is 
the total NOx flowrate produced throughout the mixer 
including any that may be produced upstream of the orifice 
injection due to the upstream jet flow recirculation in some 
of the configurations. In order to make relative comparisons 
between the mixing and NOx, the SUM-B and CUM-NOx 
parameters have been scaled to adequately fit in the bar 
charts ranging from 0 to 0.35 and all configurations are 
equally scaled to eliminate any bias. Note that these charts 
are only useful in the extrapolation of general trends and 
relative comparisons between configurations and should not 
be used for absolute evaluations. The results indicate that 
even though the round orifice has stronger jet penetration and 
lower cumulative NOx production as compared to the slanted 
slots which generally underpenetrate. Figure-7 indicates the 
location of NOx production for the round hole and elongated 
slot (L/W=8). The slot configuration (#24), shown in the 
interpolated iso-surface for max NOx production in 
Figure-7b indicates that additional NOx is produced below 
the orifice near the center of the mixer due to shallow pen- 
etration and induced swirl (see Figure-7d). Also note that the 
results generally indicate an increase in NOx production near 
the x/R=l for both best and worst NOx production configura- 
tions shown in Figure-7, which is caused by mixing the hot 
gases entering the domain and the cool gases from the jets 
along with the high residence time which mainly occurs 
downstream of the orifice. 

The mixing results indicate marginal improvement in the 
mixing nonuniformity with increased orifice aspect ratio 


(best mixer is the L/W=8 slot conf # 24) as shown consis- 
tently by all the three statistical parameters (AMIX, MMIX, . 
and volume histograms) in Figure-5a. Both mixing results 
and NOx predictions show no general correlating trends. 
Note that the volume histograms denoting the mixture 
non-uniformity indicate that the round orifice has a signifi- 
cantly large portion of the overall volume at the nondimen- 
sionalized equivalence ratio (f) near zero (see Figure-7e) as 
compared to the L/W=8 slanted slot in Figure-7f. This 
indicates that the NOx production with the round orifice is 
low since the jet flow quenches the mainstream flow. It is 
also recognized (although not analyzed in this study) that the 
production of CO and unburned hydrocarbons may be high 
since the jet flow appears to be overly quenching the rich 
zone products to a value much below the equilibrium value. 
It is also recognized that the unmixed jet could eventually 
mix with the mainstream flow downstream of the x/R=l 
boundary and therefore produce NOx. 

J=52 

The results at J=52 are shown in Figures-5b, 8, and 9. At 
J=52, the increase in penetration for the round jet appears to 
increase the NOx production since the wake behind the jet is 
larger (see top plot in Figure-8). Similarly, the slanted 
orifices appear to have better NOx performance as compared 
to the round orifice due to a more shallow jet penetration. 
Therefore, the jet penetration appears to be an essential 
parameter in controlling the NOx production. The lowest 
NOx configuration was the 45 degree slot L/W=4 (conf #29) 
and the highest NOx producer at J=52 is the L/W=8 slot. 
Both of these configurations are shown in Figure-9. 

The area weighted parameters (AMIX, MMIX, and volume 
integrals) indicate improvement in mixing uniformity as the 
slot aspect ratio increases (assume the round hole as a slot 
with an aspect ratio L/W=l). Both NOx and mixing trends 
are seemingly opposing each other in this comparison since 
the best mixing configuration is the highest NOx producer 
and these results seem to agree with the comparison at J=25 
summarized previously. The volume histograms shown in 
Figure-9e and f however indicate marginal change in the 
overall mixing flowfield equivalence ratio non-uniformity. 

J=80 

The results at J=80 are shown in Figures-5c, 10, and 11 As J 
is increased to 80, the results indicate that the highest NOx 
production is attained with the round holes and lowest NOx 
production is attained with the 45 degree slots, L/W=8. The 
ratio of the highest to the lowest NOx production is about 1.8 
(see Figure-5c). The round orifice jet appears to penetrate to 
the core of the mixer as shown in Figure-10. The 45° slot 
with L/W=4 seems to be somewhat over penetrating but the 
L/W=8 slot appears to approach optimum penetration. The 
best and worst NOx production configurations at J=80 for the 
change in orifice geometry are shown in Figure-1 1. The 45 
degree slot with aspect ratio (L/W=8) indicates no NOx 
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formation upstream of the orifice (see Figure- lie) as 
opposed to the over-penetrating jet produced by the hole 
(Figure-1 Id). Both configurations however produce NOx 
downstream of the orifice near the outer wall behind the 
orifice, it is also worthy to note that the secondary peak at 
f=0.31 shown in Figure- Ilf corresponds to significant NOx 
production for the round hole. 

The NOx production reduction trends appear to agree with 
the improvement in the mixing at J=80 as one would expect. 
This probably is due to the general decrease in jet penetration 
as the slot aspect ratio is increased from 1 (round) to 8 only 
superseded by the fact that optimum jet penetration is being 
approached from the over-penetrating side. Conversely, the 
previous comparison at J=25, optimum jet penetration 
appears to be approached from the under-penetrating side, 
and therefore NOx production and mixing trends appear to 
contradict. It is therefore apparent that optimizing the jet 
penetration could lead to minimized NOx formation as long 
as the downstream region behind the orifice is small enough 
not to produce NOx. Consequently, deviations from this 
optimum jet penetration could lead to increased NOx forma- 
tion. 

Figure- 12 indicates the comparison of the NOx production 
evolution throughout the mixing section for the change in 
orifice configuration at J=25, 52, and 80. The results show 
significant NOx production downstream of the orifice at x/R 
> l which is important to note. Several cases that 
over-penetrate indicate NOx production upstream of the 
orifice (see Figure- 12e and f). It is also important to note that 
the lowest overall NOx production configuration of the 9 
cases cited in this section is the 45 degree slot L/W=8 at J=80 
(see Table-3). 

Effect of .1 on NOx production and mixing non-uniformity 

Three different orifice shapes were compared for increasing 
J: round orifices (configuration numbers 19, 23, and 22), 45 
degree slanted slots with aspect ratio L/W=4 (configuration 
numbers 20, 29. and 30), and 45 degree slanted slots with 
aspect ratio L/W=8 (configuration numbers 24, 25, and 26) 
see Figure-2. 

Round orifices 

The effect of the increase in J on the NOx production of the 
round orifices indicates that over-penetration tends to 
increase the NOx formation by increasing the orifice wake 
volume as shown in Figures- 13a. 14, and 15. The increase in 
jet penetration produces an upstream recirculation near the 
center core of the mixer (see Figure- 14) and this added 
mixing structure produces NOx formation upstream of the 
orifice (see Figure- 13a). The area weighted parameters 
AMIX and MMIX indicate increased mixing non-uniformity 
for increasing J. The best mixing configuration is the lowest 
NOx producer in this case (Round hole conf # 19 at J=25) 
and worst mixing configuration is the round hole at J=80 due 
to the over-penetrating jet shown in Figure- 15. 


45° Slots, with aspect ratio L/W=4 

The effect of the increase in J on the NOx production of the 
45° slots with aspect ratio L/W=4 is shown in Figures-13b. 
16 and 17. The results in Figure- 12b indicate that the NOx 
production with this slanted slot decreases with increasing J 
from 25 to 80 probably mostly due to the optimization of the 
jet penetration as approached from the under-penetrating 
side. The planar mixing non-uniformity deviations (AMIX 
and MMIX) indicate deteriorating mixing performance at 
x/R=l with increasing J from 25 to 80 which generally 
opposes the NOx production trends (see Figure- 13b). This is 
caused by the hot mainstream flow which tends to be strati- 
fied near the outer wall region of the mixer. This phenomena 
was discussed previously by Oechsle, Mongia, and Holdeman 
(1993). The best mixing configuration (45 degree slot. 
L/W=4 at J=25) appears to be the highest NOx producer. 
Similarly, the worst mixing configuration (45 degree slot. 
L/W=4 at J=80) produces the lowest NOx producer. 

45° Slots, with aspect ratio L/W=8 

The effect of the increase of J on the NOx production of the 
45° slots with aspect ratio L/W=8 is shown in Figures-13c. 
18, and 19. The results indicate that this configuration 
approaches optimal jet penetration from the 
under-penetrating side with the increase in J from 25 to 80. 
Minimum NOx production is obtained with the highest J 
condition and highest NOx production obtained at J=25 (see 
Figure-I3c). Figures-19b and d show the effect of jet 
under-penetration on NOx formation near the center of the 
mixer for configuration # 24 (L/W=8 slot, J=25). As J is 
increased to J=80, the jet penetration sufficiently prevents 
any NOx production in the center core of the mixer 
(Figure- 19a and c). The planar deviations (AMIX and 
MMIX) indicate small changes in mixing non-uniformity for 
increasing J probably due to the weak effect of jet penetra- 
tion with changes in J for this particular orifice shape as 
shown in Table-3. Aside from a slightly higher peak at f=0.1, 
the volume histograms shown in Figures-19e and for the best 
and worst NOx production cases are nearly the same. This 
also correlates with the planar mimxing deviation parameters 
AMIX and MMIX results shown in Figure- 13c. 

Figure-20 shows the history of the NOx formation throughout 
the mixer for the different orifice shapes as J is increased 
from 25 to 80. Note that the rate of increase in NOx cumula- 
tive production (Figure-20b and d) are nearly the same for 
the round annd 45 degree slots (L/W=4) up to x/R=l. 

Effect of slot slant angle on NOx production and mixing 
non-uniform itv 

In this section, a parametric study was performed by varying 
the slot slant angle from 0° (in line with the flow) to 90° 
(transverse direction) for the aspect ratio L/W=4 slot only. 
The results are presented for constant J conditions. The cases 
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compared are: at J=25, configuration numbers are 34, 27, and 
20 at slant angles 0, 22.5, and 45 degrees respectively: at 
J=52, configuration numbers are 21, 28, 29, and 31 at slant 
angles of 0, 22.5, 45, and 67.5 degrees respectively: and at 
J=80, configuration numbers are 32, 35, 30, 33, and 36 for 
slant angles of 0, 22.5, 45, 67.5, and 90 degrees respectively. 
Note that some of the slant angles have been omitted from 
the study since they do not fit in the mixer. 


J=25 

At the lowest J condition (J=25), the change in NOx produc- 
tion for slant angles (varying from 0 to 45°) are shown in 
Figures-21a, 22, and 23. It is apparent that slot slant angle at 
low J has a significant effect on NOx production at x/R= 1 
with a 45% increase in NOx production from the best to 
worst (see Figure-21a). In addition, the jet penetration also 
changes significantly in these three cases (see Figure-22 and 
Table-3). The trends of both NOx formation and mixing 
appear to contradict since minimum NOx formation was 
obtained with the configuration with a the highest jet pen- 
etration which is the aligned slot. However, improved mix- 
ing was obtained with the 45 degree slot. All statistical 
mixing parameters indicate that the 45 degree slot is the most 
uniformly mixed of all three cases as shown in Figure-21a. 
The best and worst NOx production configurations for the 
variation of the slot slant angle at J=25 are plotted in 
Figure-23. The mixing results shown in Figure-23e indicate 
that the (0 degree) aligned slot has significant portion of the 
volume near f=0 as opposed to Figre-23f for the 45 degree 
slot. This means that NOx production tends to be lower for 
the slot however, this also means that mixing may eventually 
take place downstream of the mixer and NOx may be pro- 
duced. 

J=52 

The results for the comparison at J=52 are shown in 
Figures-21b, 24, and 25. At J=52, the NOx production 
appears to change significantly with slant angle and the 
lowest NOx producer was obtained with the largest slant 
angle slot (67.5°) slant angle slot as seen in Figure-21b. 
Conversely, the highest NOx (only 25% above the lowest 
NOx producer) was obtained with the 0° slant (in-line slot). 
Note that the NOx production appears to be directly propor- 
tional to the increase in slant angle, which itself is inversely 
proportional to the jet penetration (see Figure-24). The 
mixture of gases which recirculate upstream of the orifice 
tend to produce NOx which would not be produced in a 
configuration producing more optimum jet penetration as 
shown for the 67.5 degree slanted slot results (configuration 
# 31 in Figure-25c and d). 

The area and mass flow weighted mixing non-uniformity 
results shown in Figure-2 lb indicate that slot slant angle 
variation from 0 to 45 degrees does not appear to signifi- 
cantly affect the mixing performance. There however seems 


to be a noticeable improvement in the mixing when the slant 
angle is increased from 45 to 67.5 degrees. All statistical 
parameters indicate the same trend. In addition, the histo- 
gram plots in Figure-25e and f indicate that the lowest NOx 
configuration gives a histogram shape that would more likely 
indicate better equivalence ratio uniformity with a peak near 
fequjt > see Figure-25e) as opposed to the worst NOx produc- 
tion configuration identified in Figure-25f. Apparently, the 
blockage value and jet penetration are optimal such that only 
a small amount of the mainstream flow is entrained in the 
wake area behind the orifice. Perhaps the combination of the 
jet penetration and blockage forms the right size and shape 
wake region to optimize mixing near the outer wall of the 
circular mixer. 

J=80 

The results for the comparison at J=80 are shown in 
Figures-21c, 26, and 27. the NOx production changes signifi- 
cantly with slant angle and the lowest NOx was obtained 
with the 67.5° slant angle slot at x/R=l as seen in 
Figures-21c. It is also important to note that the 67.5° slot 
produced the lowest NOx compared to all the 18 configura- 
tions reported herein and this correlates with one of the best 
mixing configuration as reported by Oechsle. Mongia, and 
Holdeman (1992). As the slot slant angle increases, mixing 
improves and NOx production decreases. The highest NOx 
production configuration at J=80 for varying slant angles is 
the aligned slot (see Figure-27b, d, and f). Figure-27 shows 
that NOx formation around the orifice becomes significant 
for the low slant angle high jet penetration configuration #32 
(see Figures 27b and d). As the slant angle increases, the 
NOx production becomes localized only in the region behind 
the orifice (see Figures-27a and c). It is therefore evident 
that the jet penetration, which is directly proportional to the 
slant angle, affects the NOx production. Adequate penetra- 
tion is therefore necessary to minimize NOx and deviations 
from this value (for any given configuration) would probably 
tend to increase this emission in a mixing zone. However, it 
is also evident from the results shown in Figure-27c that 
other factors, in addition to the jet penetration, generally 
affect the NOx production and this will be examined later in 
this paper. 

The NOx production comparison in Figure-28 indicate that 
most of the configurations except for the 67.5 degree slanted 
slot (L/W=4) shown in Figure-28e and f at J=80 produce 
significant NOx downstream of the x/R=l boundary. There- 
fore, configurations with moderate penetration with good 
mixing near the outer wall of the mixer generally produce 
low NOx in a circular mixer. Over-penetrating jets produce 
NOx upstream of the orifice. However, the amount of NOx 
produced upstream of the orifice as shown in Figure-28 f is 
generally low compared to that produced in the mixing 
section itself between x/R=0 to 1. 

According to the mixing non-uniformity results, the best 
mixing configuration is the transverse slot configuration #36. 
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This configuration also is the best overall mixing configura- 
tion of all 18 compared in this paper. The trends of both the 
mixing results and NOx results do not appear to correlate in 
this comparison. Worst mixing was obtained with the 45 
degree slot at J=80 due to the flow stratification similar to 
that described by Oechsle, Mongia, and Holdeman (1992). 
Therefore it is apparent from this comparison as well as the 
previous ones described above that NOx production does not 
correlate with the statistical parameters used to describe the 
mixing criteria. Therefore selecting a specific mixer configu- 
ration for low NOx production just on the basis of the mixing 
non-uniformity flowfield requires a different approach. The 
results of the NOx and mixing criteria correlations are 
described in the next section. 


Configuration optimization at different J conditions 


The best and worst NOx production cases as a function of J 
are shown in Figure-29. AH configurations at a particular J 
condition are compared, noting that both MR and DR are 
constant in the comparison. It is also worth considering that 
the 18 configurations analyzed herein may not display the 
particularly optimum condition for either best mixing (best 
jet penetration and/or NOx production) for a particular J 
since the number of orifices is not considered as a variable in 
this study. It is apparent in Figure- 29 that the bulk of the 
NOx production is formed in the wake and downstream 
region behind the orifice and these results are also shown in 
Figure-30. In this figure, the localized NOx production for 
aH configurations is plotted as a function of non-dimensional 
position. The ±2 ct variations indicate the spread of the data 
with respect to the numerical mean as a function of position. 
Figure-25a indicates that the NOx production increases with 
distance from the orifice leading edge, Figure-25b indicates 
that NOx formation is more concentrated towards the outer 
wall of the mixer, and finally, Figure-25c indicates that the 
NOx is generally independent of transverse position in this 
global comparison. 


The round hole appears to be lowest NOx production con- 
figuration at the lowest J condition due to its optimum pen- 
etration. The optimum number of orifices for best penetra- 
tion is also calculated by Holdeman (1993) using Eq-6 where 
n is the number of orifices and C=2.5. Based on Eq-6, the 
optimum number of holes is 9 which is similar to the 8 round 
hole case modeled in conf #19. Similarly for J= 52 and 80, 
the optimum number of round holes would be 13 and 16 
respectively, which is consistent with the observed 
over-penetration for 8 holes at these momentum-flux ratios. 


jr\T 2J 

n = 

C 


( 6 ) 


For 8 orifices, increasing the operating condition (J) appears 
to shift the optimum low NOx configuration from the round 
holes to the slanted slots. Lowest NOx was produced with 
the large slant angle slots (67.5 degrees) with aspect ratio of 
4. It is believed that the optimum L/W=8 slot was not within 
the J range analyzed in this parametric study since these 


exhibited the lowest jet penetration. Therefore, if RQL 
applications require a high mixing region pressure drop (J 
above 80), it is possible that the optimum low NOx configu- 
ration could be shifted to this type of orifice. This extrapola- 
tion is based on the fact that the second lowest overall NOx 
production configuration of those reported herein is the 
L/W=8 slot at 45 degree slant angle at J=80. 


Correlation between NOx production and mixing results 

After reviewing these results, it becomes apparent that jet 
wake volume is a possible correlating parameter between 
mixing and NOx production. For this purpose, a 
non-dimensional wake volume parameter was constructed 
using Eq-7. This non dimensional volume approximates the 
volume of the wake behind the orifice up to x/R=l. Note that 
it is only an approximation because this parameter does not 
consider the jet penetration transition angle as the jet flow is 
assimilated into the mainstream flow. Therefore its purpose 
is only to establish general trends between mixing flow 
characteristics and the cumulative NOx production rate up to 
x/R=l. 


Wake Vol 


■1 


1-(1 -pen) 


(blockage) (1 - (T.E. x/R)) 


(7) 


Where, 


pen = non-dimensional jet penetration (r/R) from th 

outer wall of the mixer (see Table-3) 
blockage = orifice blockage (see Tables- 1, 2, and 3) 

T.E. x/R = orifice trailing edge x/R, (see Tables- 1, & 2). 

The results showing the correlation between NOx production 
and the non-dimensional wake volume are given in 
Figure-3 1. The results shown in Figure-3 la indicate that as 
the non-dimensional shape orientation ratio (blockage / T.E. 
x/R) increases, the NOx production decreases regardless of J 
for the symmetrical orifice configurations. The 
non-dimensoinal orifice shape orientation ratio is the small- 
est for the aligned slot and the largest for the transverse slot 
(conf # 36). The results shown in Figure-3 lb indicate that as 
slot slant angle increases, NOx increases only for the low J 
condition since a somewhat under penetrated case (22.5 
degree slot conf # 27) becomes even more under penetrated 
with the increase of slant angle to 45 degree (conf # 20). 
This trend however reverses for the higher J conditions which 
have sufficient jet momentum to allow for slot slant angle 
variation to a more optimized geometric configuration. All 
the data points plotted in Figure-3 lb are only for the L/W=4 
slots. The results shown in Figure-3 lc indicate the approach 
towards NOx reduction optimization of each orifice configu- 
ration as J is increased. For the orifices with lower orifice 
shape orientation ratio such as the aligned slot, 22.5 degree 
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slot, and the round hole, the NOx production increases for 
increasing J, generally as a result of over penetration. How- 
ever, higher orifice shape orientation ratios such as the 67.5 
degree slot and transverse slot reverse the trend due to their 
inherent low jet penetration, and approach the optimum 
desired jet penetration from the under-penetrated side as J 
increases from 25 through 80. This correlation indicates that 
the non-dimensional wake volume is a possible correlating 
parameter linking the mixing flow structure to the NOx 
production; however, improvements could be made if param- 
eters such as J, orifice shape orientation, and velocity distri- 
bution are added in order to obtain a more generalized corre- 
lating parameter. 

Conclusions 

1) The mixing statistical parameters in most cases do not 
correlate with the NOx production rates at x/R=l. The planar 
variances at x/R=l lack the mixing history throughout the 
entire mixing region and the volume integrals require the 
residence time factor for completeness. 

2) NOx production is shown to be highly related to the jet 
penetration. Over penetrating configurations show increased 
NOx production as so do under-penetrating cases. At low J 
conditions, an orifice with large orientation shape ratio tends 
to under penetrate and show characteristics of a non-optimal 
NOx reduction configuration: and subsequently higher J, 
tends to decrease the NOx production. Similarly, in the 
range of J analyzed herein, orifices with low orientation 
shape ratios tend to over penetrate indicating non-optimal 
mixer configurations; and subsequent increase J, increases 
the amount of NOx produced. 

3) NOx production correlation with position indicates that for 
all the 18 configurations summarized in this report, most of 
the NOx produced occurs near or in the orifice wake region. 
After correlating the NOx production with the 
non-dimensional wake volume, results indicate that trends do 
exist but additional terms such as orifice shape ratio, veloc- 
ity, and J need to be incorporated to obtain a more general 
relation to correlate mixing with NOx. 

4) At the lowest J condition, the circular orifice indicated the 
lowest NOx production, and at higher J conditions, the 
optimization shifted towards the large slant angle slots such 
as the 67.5 degree slot L/W=4. High NOx production was 
obtained with the aligned slots due to the over-penetrating jet 
structure. Over penetrating jets also cause additional NOx to 
be produced upstream of the orifice. 

5) The lowest NOx formation in a mixing region may not 
necessarily mean lowest NOx in the overall RQL combustor. 
The mixture non-uniformity in the mixing section needs to 
also be evaluated along with the NOx production to prevent 
any further unmixed flow from developing in the lean zone 
and producing NOx where is currently not accounted for in 
this paper. 
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PARAMETER 

UNITS 

19 

20 

21 

22 

23 

24 

25 

26 

27 

OVERALL 











J (jet/mian) 


26.7 

30.5 

51.1 

84.1 

553 


50.9 

88.4 

28.0 

MR (jet/main) 


2.96 

2.96 

2.96 

2.96 

2.96 

2.96 

2.96 

2.96 

2.96 

DR(jet/main) 


2.278 

2378 

2.278 

HSSSii 

221% 

2.278 

221% 

2378 

2.278 

VR (iet/main) 


3.424 

3.659 

4.736 

6.080 


3312 


6333 

3.506 

Ar (iet/main) 


0380 

0.355 

0214 

0.214 

0.264 

0.370 



0.371 

Phi (rich zone) 


L80 ~ ' 

1.80 


1.80 


1.80 

1.80 


130 

Phi (lean zone) 


0.42 

0.42 


0.42 

HEZ9I 

0.42 


0.42 

0.42 

DP/P 

% 

138 

1.58 

mum 

435 

2.86 

1.45 

'WKfwm 

4.57 

1.45 












MAINSTREAM 










MHIM 

P main 

atm. 

14.1 

14.1 

14.1 

14.1 

14.1 




KIM 

T main 

K 


mem 

HI 

■Mil 

mmm 



WBSM 

KUjH 

U main 

m/sec 

25.4 


25.4 

25.4 

25.4 

| 

25.4 

25.4 


m main 

mamm 

0.266 

0266 

0.266 

0366 

0366 

| 

0366 

0366 


Mixer diameter 


0.076 

0.076 

0.076 

0.076 

0.076 


0.076 

0.076 













JET 


ROUND 

SLOT 

SLOT 

ROUND 

ROUND 

SLOT 

SLOT 

SLOT 

SLOT 

P let 

atm. 

14.1 

14.1 

14.1 

14.1 

14.1 

14.1 

14.1 

14.1 

14.1 

_ — T. J:- 

T iet 

K 

950 

950 

950 

950 

950 

950 

950 

950 

950 

V iet 

m/sec 

86.8 

92.8 

120.1 

154.1 

124.1 

89.0 

119.8 

158.0 

88.9 

J 

m iet 

kg/sec 

0.788 

0.788 

0.788 

0.788 

0.788 

0.788 

0.788 

0.788 

0.788 

Jtz 

ACd iet/row 

m 2 

1.732e-03 

1.620e-03 

1352e-03 

9.753e-04 

1.202e-03 

1.688e-03 

1.254e-03 

9.512e-04 

1.691e-03 

Orifice Cd 


0.75 

0.74 

0.77 

0.74 

0.76 

0.75 

0.77 

0.76 

0.75 

Number of orifices 


8 

8 

8 

8 

8 

8 

8 

8 

8 

Orifice L/W 


1 

4 

4 

1 

1 

8 

8 

8 

4 

Orifice slant angle 


0 

45 

0 

0 

0 

45 

45 

45 

22.5 

Orifice length (L) 

m 

0.017 

0.029 

0.026 

0.012 

0.014 

0.042 

0.036 

0.031 

0.030 

Orifice width (W) 

m 

0.017 

0.007 

0.006 

0.012 

0.014 

0.005 

0.004 

0.004 

0.007 

Blockage 


0.5548 

0.7630 

0.2149 

0.4164 

0.4623 

1.0351 

0.8922 

0.7770 

0.5364 

Orifice T.E. x/R 


0.4357 

0.5992 

0.6750 

0.3270 

03631 

0.8130 

0.7007 

0.6103 

0.7398 


Table-1. Configuration description and parameter specifications for configurations 19 through 27 





















































































PARAMETER 

UNITS 

28 

29 

30 

31 

32 

33 

34 

35 

36 

OVERALL 






— 


— 



J (jet/mian) 


53-2 

57.6 

92.9 

59.9 

93.0 

992 

HI 

mmm 

106.2 

MR (iet/main) 


2.96 

2.96 

2.96 

msm 

■H 

2.96 

2.96 

2.96 

2.96 

DR(jet/main) 


2.278 

2.278 

SI 

2.278 

2278 

2.278 


2278 

2.278 

VR (iet/main) 


4.832 

5.033 


5.128 

6293 



6.219 

6.831 

Ar (iet/main) 


0.269 

0.258 

■ mm 

0253 

0203 


0.385 

0209 

0.190 

Phi (rich zone) 




1.80 




120 

1.80 

1.80 

Phi (lean zone) 



■S3SI 

0.42 

mam 

HjSEH 

WESk H 

0.42 

0.42 

0.42 

DP/P ! 

% 


2.98 

4.81 

' 

3.10 

4.81 

5.13 

124 

425 

5.49 












MAINSTREAM 











P main 

— 

14.1 

WBSM 

14.1 

14.1 


14.1 


Bill 


T main 

K 

2,164 

mam 

IM 

2,164 

M&ZM 

nopi 

MMM 

KfiH 


U main 

m/sec 

25.4 


HESS 


MSSM 

HEUK 

25.4 

IKS 

25.4 

m main 

kg/sec 

0.266 

0-266 

0-266 


0266 

0266 

0.266 

0266 

0.266 

Mixer diameter 

m 

0.076 

0.076 

0.076 I 

0.076 

0.076 

0.076 

0.076 

0.076 

0.076 












JET 


SLOT 

SLOT 

SLOT 

SLOT 

SLOT 

SLOT 

SLOT 

SLOT 

SLOT 

P iet 

atm. 

14.1 

14.1 

14.1 

14.1 

14.1 

14.1 

14.1 

14.1 

14.1 

- cl 

T jet 

K 

950 

950 

950 

950 

950 

950 

950 

950 

950 

Viet 

m/sec 

122.5 

127.6 

162.0 

130.0 

162.1 

167.4 

85.6 

157.7 

173.2 

m jet 

kg/sec 

0.788 

0.788 

0.788 

0.788 

0.788 

0.788 

0.788 

0.788 

0.788 

ACd iet/row 

m 2 

1227e-03 

1.178e-03 

9.277e-04 

1.156e-03 

9274e-04 

8.977e-04 

1.755e-03 

9231e-04 

8.679e-04 

Orifice Cd 


0.75 

0.71 

0.73 

0.70 

0.74 

0.71 

0.79 

0.76 

0.68 

Number of orifices 


8 

8 

8 

8 

8 

8 

8 

8 

8 

Orifice L/W 


4 

4 

4 

4 

4 

4 

4 

4 

4 

Orifice slant angle 


22.5 

45 

45 

67.5 

0 

672 

0 

222 

90 

Orifice length (L) 

m 

0.025 

0.025 

0.022 

0.025 

0.022 

0.022 

0.030 

0.022 

0.021 

Orifice width (W) 

m 

0.006 

0.006 

0.006 

0.006 

0.006 

0.005 

0.008 

0.006 

0.005 

Blockage 


0.4569 

0.6506 

02774 

0.7787 

0.1849 

0.6863 

02544 

0.4027 

0.7156 

Orifice T.E. x/R 


0.6300 

0.5110 

0.4535 

0.3483 

0.5810 

02070 

0.7991 

02554 

0.1405 


Table-2. Configuration description and parameter specifications for configurations 28 through 36 


















































































1 

ppaEiiJiTsliiKiimil 

Conf 

AMIX 

MMIX 





T.E.x/R 

r/R 

# 

(8) x/R=l 

@ x/R=l 



mm 

i 


from outer wal 

mm 


0.2010 

0.398 

0.0113 

26.7 

0.555 

0.436 

0.425 

20 


0.1590 

0.331 

0.018 

30.5 

0.763 

0.599 

0.419 

mm 



0.418 

0.0172 

IlfOi 

■ flMKWi 

0.675 


mm 



■mg 1 

0.0186 

1 

WS-M 

msmm 

0.6 

23 


0.2450 


0.015 

ESI 



0.536 

mm 


0.1610 

0.34 

0.0179 

eh 


0.813 

WmOKiE^m 

mm 



0.352 


Kin 


0.700 


26 

0.2090 

0.1690 

Bili2 


V«T!VH 


0.610 

I 

27 

0.2340 

0.2000 


0.0123 

28 


0.740 

0.43 

28 

0.2420 

0.2000 

0.442 

0.0129 

53.2 

■SMB 

0.630 

0.517 

29 

0.2330 


0.419 

0.0129 

57.6 

USSU 

0.511 

0.536 

Wfl 


BETf 

0.475 

0.0116 

92.9 

■ 

0.454 

0.608 

31 

0.1500 

0.1690 

0.319 

■n 

59.9 




32 

0.1680 


0.34 


93 


0.581 


33 

0.2200 

0.1830 

0.403 

0.0086 

99.3 

0.686 

0.307 

0.541 

34 

0.2440 

0.2330 

0.477 

0.0125 

26 

0.254 

0.799 

0.53 

35 

0.2160 

0.1940 

0.41 

0.0155 

88.1 

0.403 

0.555 

0.587 

36 

0.1340 

0.1230 

0.257 

0.0136 

106 

0.716 

0.141 

0.557 


Table-3. Mixing non-uniformity and NOx production results for configuraions # 19 through 36. 








































































Orifice 
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Cylindrical 

domain 
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if 


x/R=- 1 4 | FLO w DIRECTION x/r= 6 


Tangential direction 


x=0 


T 

2R=76 mm 


Figure- 1. a) Geometric configuration of the mixing section, and 

b) orthogonal view of the computational grid with a 45° slot 
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Configuration # 
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90.0" 
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67.5° 
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Nominal jet-to-mainstream momentum flux ratio (J) 


Figure-2. Matrix depiction of the analyzed configurations in a circular mixer 


a. Example plot of the equivalence raito and velocity distributions 


Flame im Ike axial-radial direction through the center of the orifice 


NORMAIJZED 

♦ 

MAINSTREAM • 14 


b. Example plot of the NO x production evolution in the mixing section 


r/K 



d. Example plot of the sub-volume division in the mixer domain 


for the NO x analysis 




gr NO /l sec m*) 


Edge closest lo 



c. Example plot of the iso-contour of constant NO t production at 10 3 
gr of N0 X / {sec m 3 /. Mixer domain is 0 < xJR < I 


Edge closest to 
the observer 


Figure-3 . Example plots for the observational analysis 



Figure-4. Description of the statistical analysis for the mixing 
non-uniformity evaluation 
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Figure-5. Effect of the change in orifice shape on the 
mixing non-uniformity and NOx production 
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Equivalence ratio 


velocity 


Figure-6. Effect of the change in orifice geometry on the equivalence ratio and velocity distributions 
Top to bottom: round , 45° slot with L/W-4, and 45 0 slot with LfW=8 
7=25, MR=2.96, DR- 2.28, 8 orifices/row 
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BEST CONFIG URA TION WORST CONFIG URA TION 

(a, b) Interpolated N0 X production isopleth at I0 3 gr of NO J {sec m 3 ) 

(c, d) NO x production development through the mixing section 

(e,f) volume histograms depicting the mixing non-uniformity for the mixing domain (0 < xJR < l) 

Figure-7. Comparison of the best ( left column) and worst (right column) NO t production configurations 
with variation of the orifice configuration at J-25. 
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Figure-8. Effect of the change in orifice geometry on the equivalence ratio and velocity distributions 
Top to bottom: round , 45° slot with L/W=4> and 45 0 slot with L/W=8 
J=52 y MR=2.96, D/?=2.2S, 8 orifices/row 
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BEST CONFIGURATION WORST CONFIGURATION 

(a, b) Interpolated NO x production isopleth at I0 3 gr of N0 X / { sec m 3 } 

(c, d) NO x production development through the mixing section 

(e,J) volume histograms depicting the mixing non-uniformity for the domain (0 < xJR < 1) 

Figure-9. Comparison of the best ( left column) and worst (right column) NO x production configurations 
with variation of the orifice configuration at J-52. 
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Figure-] 0. Effect of the change in orifice geometry on the equivalence ratio and velocity distributions 
Top to bottom : round , 45° slot with L/W=4 , and 45° slot with L/W=8 
J=80 , MR- 2.96, DR=2.28 , <S orifices/row 
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45 1 slot, UW= 8, J= 88.4, MR= 2.96, DR= 2.28 b. 


Round hole, J=84.1, MR=2.96. DR=2.28 



BEST CONFIG U RATION WORST CONFIG URA TION 

(a, b) Interpolated NO x production isopleth at lit' gr of NO x / {sec m’j 
(c y d) NO x production development through the mixing section 

(tyj) volume histograms depicting the mixing non-uniformity for the domain (0 < x/R < 1) 

Figure-11 . Comparison of the best (left column) and worst (right column) NO production configurations 
with variation of orifice configuration at 7 = 80.. 
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a. 



x/R DISTANCE FROM THE LEADING EDGE OF THE ORIFICE 


b. 



c. d. 




x/R DISTANCE FROM THE LEADING EDGE OF THE ORIFICE x/R DISTANCE FROM THE LEADING EDGE OF THE ORIFICE 



xJR DISTANCE FROM THE LEADING EDGE OF THE ORIFICE 



(a, b) Effect of the change in orifice shape on NO x production at J= 25 
(c, d) Effect of the change in orifice shape on NO production at J=52 
(e,f) Effect of the change in orifice shape on NO production at J=80 

Figure-12. NO x production development throughout the mixing section. NO x planar production (left column) 
and NO x cumulative production ( right column ). 
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Figure- 13. Effect of the change in J on the mixing non-uniformity 
and NOx production rate. 
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Figure-8. Effect of the change in orifice geometry on the equivalence ratio and velocity distributions 
Top to bottom: round , 45° slot with L/W=4> and 45 0 slot with L/W=8 
J=52 y MR=2.96, D/?=2.2S, 8 orifices/row 
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BEST CONFIG URA TION WORST CONFIG URA TION 

(a, b) Interpolated NO x production isopleth at 10 3 gr of NO x / (sec m 3 } 

(c, d) NO x production development through the mixing section 

(e,f) volume histograms depicting the mixing non-uniformity for the domain (0 < xJR < I) 

Figure-15. Comparison of the best (left column) and worst (right column) NO x production configurations 
with variation of J for the round holes. 
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Figure-16. Effect of the change in J on the equivalence ratio and velocity distributions 

Top to bottom: 7= 25, 52, an*/ 50 
45* s/ote, L/W=4, MR=2.96, DR=2.28 , 5 orifices/row 
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BEST CONFIG ORATION WORST CONFIG URA TION 

(a, b) Interpolated NO x production isopleth at l(f gr of N0 X / {sec m 3 ) 

(c, d) N0 X production development through the mixing section 

(e,f) volume histograms depicting the mixing non-uniformity for the domain (0 < xJR < 1) 

Figure-17. Comparison of the best (left column) and worst (right column) NO x production configurations 
with variation of J from 25 to 80 for the 45 * slots with an aspect ratio L/W= 4. 
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Figure- 18. Effect of the change in J on the equivalence ratio and velocity distributions 

Top to bottom: 7= 25, 52, and 80 
45* slots, L/W=8, MR= 2.96, DR=2.28, 8 orifices/row 


45 


4? slot, LfW=8, J=88.4, MR- 2.96, DR=2.28 b. 45° slot, LTW=8, J=28.1, MR=2.96, DR=2.28 



BEST C ON FIG URA TION WORST CONFIG URA TIO N 

(a, b) Interpolated NO, production isopleth at I0 3 gr ofNO x / (sec m 3 } 

(c, d) NO gl production development through the mixing section 

(e,f) volume histograms depicting the mixing non-uniformity for the domain (0 < xJR < l) 

Figure-19. Comparison of the best (left column) and worst (right column) NO x production configurations 
with variation of J from 25 to 80 for the 45* slanted slots with aspect ratio of L/W—8 
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a. 



xJR DISTANCE FROM THE LEADING EDGE OF THE ORIFICE 


c. 



e. 



b . 



xJR DISTANCE FROM THE LEADING EDGE OF THE ORIFICE 


d. 



x/R DISTANCE FROM THE LEADING EDGE OF THE ORIFICE 



xJR DISTANCE FROM THE LEADING EDGE OF THE ORIFICE 


(a,b) Effect of the change in J for the round holes 

(c y d) Effect of the change in J for the 45 } slots with aspect ratio (L/W) - 4 
(e y f) Effect of the change in J for the 45 } slots with aspect ratio (L/W) = 8 

Figure-20 . NO x production development throughout the mixing section . NO x planar production ( left column ) 
and NO x cumulative production ( right column ). 
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Figure-21. Effect of the change in slot slant angle on the mixing 

non-uniformity and NOx production rate for the L/W=4 slot 
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Figure-22. Effect of the change in slot slant angle on the equivalence ratio and velocity distributions 
Top to bottom: slant angles 0° (in-line), 22.5°, and 45° 

7=25, slot aspect ratio L/W=4 , MR- 2.96, DR- 2.28, 8 orifices/row 
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a. 0* (aligned slot), UW=4, J=26.0, MR=2.96, DR=2.28 b. 


45* slot, UW-4, J-30.5, MR=2.96, DR=2.28 



BEST CONFIG URA TION WORST CONFIG URA TION 

(a, b) Interpolated NO x production isopleth at 10 3 gr of NO x / { sec m} 

(c, d) NO t production development through the mixing section 

(e,J) volume histograms depicting the mixing non-uniformity for the domain (0 < xJR < 1) 

Figure-23. Comparison of the best (left column) and worst (right column) NO x production configurations 
with variation of slot slant angle from o* to 45° for the aspect ratio slot of L/W= 4, at J= 25 
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Figure-24. Effect of the change in slot slant angle on the equivalence ratio and velocity distributions 
Top to bottom: slant angle- (f (in-line), 22.5*, 45*, and 67.5° 

J=52, slot aspect ratio L/W-4, MR-2.96, DR- 2.28, 8 orifices/row 
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67.5 * slot, L/W=4, J-59.9, MR=2.96, DR=2.28 b. 0 9 (aligned slot), UW=4, J=51.1, MR=2.96, DR=2.28 



BEST CONFIGURATION WORST CONFIGURATION 

(a, b ) Interpolated N0 J production isopleth at I & gr of NO x / { sec m 3 } 

(c, d) NO x production development through the mixing section 

(e,f) volume histograms depicting the mixing non-uniformity for the domain (0 < xJR < I) 

Figure-25. Comparison of the best ( left column) and worst (right column) NO x production configurations 
with variation of slot slant angle from o° to 67.5 6 for the aspect ratio slot of UW-4, at J=52. 
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Figure-26. Effect of the change in slot slant angle on the equivalence ratio and velocity distributions 
Top to bottom: slant angles 0 * (in-line), 22.5*, 45*, 67.5*, and 9<f (transverse to the flow) 
J=80, slot aspect ratio L/W-4, MR-2.96, DR=2.28, 8 orifices/row 
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BEST CONFIG URA TION WORST CONFIG URA TION 

(a, b) Interpolated NO M production isopleth at I O 1 gr of NO i / (sec m 3 } 

(c, d) NO M production development through the mixing section 

(e,f) volume histograms depicting the mixing non-uniformity for the domain (0 < x/R < I) 

Figure-27. Comparison of the best (left column) and worst (right column) NO production configurations 
with variation of slot slant angle from o* to 90° for the aspect ratio slot of L/W=4, at J=80. 
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(a, b) Effect of the change in slot slant angle on NO ^formation at J=25 
(c, d) Effect of the change in slot slant angle on NO x , formation at J=52 
(e,f) Effect of the change in slot slant angle on NO x formation at J=80 

Figure-28. NO x production development throughout the mixing section. NO x planar production (left column ) 
and NO x cumulative production ( right column). 
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a. Round hoU, 2=26.7. MR=2.96. 1)R=2.28 



e. 67 j" slot. 1JW=4, 2=59.9. MR=2.96. DR=2.28 



Z 



z 



t. 67 y slot, UW=4, J=99.3, MR=2.96, DR=2.28 



b. 45* slot, L/W=8. 2= 28. 1, MR=2.96, DR= 2.28 



i. Aligned slot, UW=4, 2=51.1, MR=2.96. DR=2.28 



f. Aligned slot, UW=4, 2=93.0, MR=2.96, DR=2.28 


Figure-29. Comparison of the best ( left column) and worst (right) local NO production 
configurations for nominal J=25 (a, b), 52 (c, d), and 80 (e,f). The interpolated surface shown 
is the 10 3 gr NO t /{sec m 3 } value. Domain is 0 < xJR < 1. 
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Figure-30. Correlation of the localized NO x production in (gr ofNO x / { sec m }) 

with position for all configurations # 19 through 36. Mixer domain is from 
0 < x/R < 1. 
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N0 X vs geometry correlation for different orifice shapes 
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Figure-31. Correlation between NO production and orifice geometry for 
configurations (#18 through 36) 
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